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Abstract: In this paper we present the evolution of a low mass model (initial mass 
M=1.5 Mq) with a very low metal content (Z = 5 x 10"^, equivalent to [Fe/H]= -2.44). 
We find that, at the beginning of the AGB phase, protons are ingested from the envelope 
in the underlying convective shell generated by the first fully developed thermal pulse. 
This peculiar phase is followed by a deep third dredge up episode, which carries to the 
surface the freshly synthesized ^'^C, ^^N and ^Li. A standard TP-AGB evolution, then, 
follows. During the proton ingestion phase, a very high neutron density is attained and 
the s-process is efficiently activated. We therefore adopt a nuclear network of about 700 
isotopes, linked by more than 1200 reactions, and we couple it with the physical evolution 
of the model. We discuss in detail the evolution of the surface chemical composition, 
starting from the proton ingestion up to the end of the TP-AGB phase. 
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1 Introduction 



The interpretation of the already available and still 
growing amount of data concerning very metal-poor 
stars fsee lBeers fc Christlieb I l|2005l ) for the classifica- 
tion of these objects) require reliable stellar models. A 
significant fraction (~25%) of stars with [Fe/H]<-2.5 
show carbon-enriched atmospheres (Carbon Enhanced 
Metal Poor stars - CEMP stars). Moreover, many 
CEMP stars (~75%) are characterized by s-process 
enhanced patterns (CEMP-s): in this paper we will 
only address to this sub-class of low metallicity stars. 
CEMP-s stars belong to the Galactic Halo and, there- 
fore, are old (~ 14 Gyr) low mass (M < 0.9 Mq) 
objects; notwithstanding, they show the concomitant 
enhancements of C, N and s-process elements, which 
are commonly ascribed to the complex interplay of nu- 
cleosynthesis, mixing and mass loss taking place along 
the Asymptotic Giant Branch (AGB). When these low 
mass stars reach the AGB, their envelope is so small 
that third dredge up (TDU) may not take place, so 
they are probably dwarfs or giants belonging to binary 
systems. Therefore, the C-enhancement and the prod- 
ucts of the neutron capture nucleosynthesis could re- 
sult from an ancient accretion process (via stellar wind 
or, less likely, via Roche lobe overflo w) from the (now 
extinct) AGB companion (se e, e.g., iRoederer et al. I 
l2008l : [Thompson et al. 1 l2008t l. Note that a recent de- 
termination of radial velocities for a sample of CEMP-s 
stars suggests t hat all of these objects are members of 
binary systems (jLucatello et al. 1l2005l ). Some of these 
stars present exotic chemical distributions and, for this 
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reason, represent the current frontier of AGB model- 
ing. 

At very low metallicities, there is a minimum mass 
(which increases with decreasing the metallicity) under 
which models suffer proton ingestion from the enve- 
lope down to the underlying convective He-intershell. 
This occurs during the first fully developed thermal 
pulse (TP). We refer to this episode as the Proton In- 
gestion Episode (PIE). Larger mass models follow a 
standard TP-AGB evolution. This feature has a l ready 
been found by many authors jHoUowelLglL-alJ [l990|; 
iFuiimoto et al . l2000l: I Su da et al. '2004':'lwamoto et all 
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2004; Stranic ro et al. 120 04: Cristallo ct al. 2007; Ca mpbeU fc Lattanzi^ 
t200a : ILau et al. 1120091 '): the occurrence of PIE is there- 
fore a robust prediction, in spite of the different physics 
adopted in the various work s. Moreover, recent 3D 
hydrodynamical simulations ([Woodward et ah"! 120081 ') 
confirm this peculiarity, which characterizes very metal- 
poor models. A full understanding of the involved 
physics is therefore mandatory. 

In next Sections, we present and discuss the evolu- 
tion of an AGB model with initial mass M=1.5 Mq 
and Z — 5 X , with emphasis on the effects on 
the stellar structure caused by the PIE. The model 
has been com puted with the FR ANEC stellar evolu- 
tionary code (jChieffi et al. I IT998I): details of the mix- 
ing scheme (see next Section) and latest improvement s 
to the code are described in ICristallo etaTI ||2009D 
We a ssume a solar-scaled initial distribution ijLodders I 
120031 ') ■ although at low metallicities a clear oxygen en- 
han cement has been obser ved by many authors (see, 
e.g., lAbia fc Reboio1ll989l ). 
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Figure 1: Evolution of the temperature within the Hc-intcrshell during the PIE. We plot the profiles at 
the beginning of proton ingestion (panel a, At=0 yr), soon before the convective shell splitting (panel b, 
At=1.6457 yr), soon after the convective shell splitting (panel c, At=1.6468 yr) and during the relaxing 
phase (panel d, At=2.1843 yr). Note how the convective shell (shaded region) in the top left panel is not 
yet fully developed. Different shading directions (panel c and d) indicate separate convective regions. See 
text for details. 



2 The Proton Ingestion Episode 



As outlined in the previous Section, at the beginning 
of the thermally pulsing phase of low mass metal-poor 
models, the convective shell powered by 3a reactions 
may extend up to the base of the H-rich envelope. 
This only occurs at low metallicities and is due to the 
flattening of the entropy peak characterizing the H- 
burning shell. At low Z, in order to compensate the 
reduction of CNO catalysts, the temperature of the H- 
shell increases and, consequently, the entropy within 
the H-shell decreases (entropy variation is defined as 
dS = dQ/T, where dQ is the exchanged heat and T 
is the temperature). Therefore, while at larger metal- 
licities mixing is pre vented by the entropy peak of the 
H-shell ijiben |[T977t ). at low metallicities protons can 
be mixed within the convective shell and captured by 
the abundant ^^C. The temperature profile within the 
convective shell runs from about T = 5.0 x 10^ K at 
the outer border to T = 2.3 x 10* K at the bottom (see 
top left panel of Fig. [!}. At such high temperatures 
H-burning is very efficient and, therefore, protons are 



captured on-flight by CNO isotopes before reaching 
the bottom of the convective shell. In these condi- 
tions, we can distinguish three different regions within 
the flash-driven convective shell: 

• the external zone, where the CNO-burning time 
scale (rcNo) is largely greater than the mixing 
time scale; 

• the inner zone, where temperature is so high 
that the mixing efficiency is largely smaller than 
the burning one; 

• an intermediate zone, where the nuclear and mix- 
ing time scales are comparable. 

The best approach to this situation is the simultane- 
ous solution of the differential equations describing the 
physical structure of a star, the nuclear burning and 
the mixing. However, it has to be noticed that this 
computational scheme can be used only with small nu- 
clear network, while it can not be applied when using 
very large network, including several hundreds of iso- 
topes and nuclear processes. For this reason, in our 
simulations we adopt a different approach. We fully 
couple the structure equations with those describing 
the nuclear burning, so that we directly take into ac- 
count the energetic feedback of nuclear reactions on the 
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Figure 2: Evolution of selected key isotopes within the He-intershell during the PIE. 



physical properties of the model (mainly the temper- 
ature profile). Moreover, this allows us to determine 
accurately the temperature gradient and, thus, the lo- 
cation and the mass extension of the zones unstable 
for convection. Hence, convection-induced mixing is 
computed b y adopting the t ime d ependent formalism 
described in IChieffl et al."1 l|200lD . where the mixing 
efhciency depends linearly on the temporal time step 
At of the model. In our simulation, during the PIE, we 
limit At to 50% of the mixing turnover time scale for 
the whole convective shell {rmix), in order to avoid the 
unrealistic fully homogenization of that region. All the 
chemical species are mixed within the whole unstable 
zone, with the exception of protons, which are mixed 
from the outer border of the convective shell down to 
the mass coordinate Mp, where tcno = /At, / < 1 
being a free parameter. According to our choice on 
At, at Mp TCNO results definitively smaller than Tmix 
(at least by a factor of 2//): this implies that protons 
consumption via nuclear burning is much more rapid 
than the supply via the convection-driven mixing. As 
a consequence, in nature below this point it is hardly 
unlikely for protons to survive. In the computation 
of PIE, we fix / = 1/3 (this implies that at Mp the 
temperature roughly corresponds to T ~ 1.3 x 10*). 
Such a choice is based on a a series of evolutionary 
sequences computed with different values of /. Our 
results show that for / > 1/3 protons are accumulated 
in cool zones where H-burning is inefficient, while for 
/ < 1/3 protons are mixed down to zones inconsistent 
with the on-flight burning. 

In panel a of Fig. [2] we report the abundances of H, 
""^^C, ^^C and at the beginning of proton ingestion. 



Note the large C abundance, due to the partial He- 
burning inside the convective shell via 3q: reactions, 
and the relatively low abundances of ^^C and ^*N, 
which result from the dilution, within the convective 
shell, of the ashes left by the H-burning shell during 
the previous interpulse phase. Similarly, in panel a of 
Fig. [3] we plot the contributions to energetics gener- 
ated by proton captures (blue long-dashed curve) and 
by Q captures (dark solid curve). Two peaks clearly 
emerge from the figure; the most internal one is due 
to 3a reactions (which power convection within the 
He-intershell) while the second peak (more external 
and shallow) comes from the CNO cycle in the H- 
burning shell. We plot apart the contribution of the 
^^C(a,n)^®0 reactions, in which we also consider the 
energy release by the related neutron captures (short- 
dashed red curve). 

Upper right panels (panels b) of Figures [TJ [5] and [3] 
depict the situation 1.6457 yr after the beginning of 
protons ingestion. Protons are captured by the abun- 
dant ^■^C, leading to the nucleosynthesis of ^^C and 
^^N. Due to the concomitant actions of mixing and 
burning, the resulting abundances strongly differ with 
respect to the ones of the radiative H-burning shell, 
in which a full CNO equilibrium is attained and the 
most abundant synthesized isotope is ^''N. During the 
PIE, the H-burning is instead always out of equilib- 
rium and this translates in a large production of ^"^C 
and, at a lower level, of ^*N, while ^^O remains prac- 
tically untouched. Technically, this is due to the fact 
that isotopes cannot locally pile up (and then reach 
the equilibrium value) due to the presence of convec- 
tive motions. Note that ^^C and ^*N are correctly 
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Figure 3: Evolution of different contributions to energetics within the He-intershell during the PIE. 



mixed down to the base of the convective shell, be- 
cause their burning timescales are always lower than 
the model temporal step. The hydrogen profile extends 
down to M ~ 0.548 M0. The very small H abundance 
under this limit basically derives from ^''N(n,p)^'*C re- 
actions: neutrons mainly come from the "C(a,n)i^O 
reaction, which is very efficiently activated at the bot- 
tom of the convective shell. This is clearly shown in 
panel b of Fig. [S] in which emerges that ^''C(a,n)^^0 
reactions (and related neutron captures) represent the 
major energy source at the bottom of the convective 
shell. At the same time, the contribution of CNO burn- 
ing boosts a second energy peak, whose mass coordi- 
nate is determined by the maximum penetration of 
protons within the convective shell. When the maxi- 
mum CNO energy peak slightly exceeds the value of 
the He-burning peak, an inversion in the temperature 
profile occurs (see panel c of Fig. [T]) and the convective 
shell splits into two sub-shells: the lower one boosted 
by the ^''C(Q,n)^®0 reactions and the upper one by the 
CNO cycle. After the splitting event, the two convec- 
tive shells follow separate evolutions. 
In the lower shell, the ^^C(Q!,n)^''0 reaction consumes 
the available ^''C and produces a very large neutron 
flux. At the bottom of the shell, a neutron density as 
large cLS Tin ~ 10^^ cm ^ is attained, allowing the syn- 
thesis of isotopes by-passed by a standard s-process. 
Since so high a neutron density moves the nucleosyn- 
thetic path away from the /3 stability valley, we ex- 
tended the nuclear network with respect to that one 
used in previous works: the number of isotopes is in- 
creased from 500 to about 700 and the number of nu- 



clear reactions grows from 700 to about 1200. 
In the upper shell, the ^^C and the ^*N continue to 
increase thanks to the burning of freshly ingested pro- 
tons from the external layers. In this phase, the CNO 
cycle is the major energy source in the star (see panel c 
of Fig. [3]l . Then, the structure starts expanding, caus- 
ing the temperature of the whole region to rapidly de- 
crease (see panel d of Fig. [1} . During this phase, both 
the energy sources boosting the two shells quench off 
(see panel d of Fig. in the meanwhile, the upper 
shell extends in mass. 

Later on, the envelope completely engulfs the upper 
shell, dredging up to the surface a large amount of ^''C 
and ^*N. The lower shell, thanks to the steep pressure 
and density gradients formed in the splitting region, 
remains part of the He-intershell region and the iso- 
topes synthesized in this region are later diluted in the 
following TP. Owing to the increase of the CN abun- 
dance in the envelope, the further evolution of this 
models resembles that of a standard AGB, preventing 
the occurrence of further PIEs. In the upper panel of 
Fig.|4]we report the temporal evolution of the mass co- 
ordinates of (top to bottom): the inner border of the 
convective envelope, the maximum energy production 
of the H-burning shell and the maximum energy pro- 
duction within the H-depleted core. In the lower panel 
we plot the H-burning and He-burning contributions to 
luminosity. Note the H-burning burst (Lh > W^Lq) 
caused by the PIE. This model undergoes 25 additional 
TDU episodes, each one followed by the formation of 
a ^^C pocket (ISt raniero et al. I l2006l : ICristaUo et al. I 
[2002), which radiatively burns during the interpulse 
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Figure 4: Upper panel: temporal evolution of the mass coordinates of (top to bottom) the inner border 
of the convective envelope, the maximum energy production of the H-burning shell and the maximum 
energy production within the H-depleted core. Note the occurrence of a deep TDU (hereafter d-TDU) 
at the beginning of the TP-AGB phase, as a consequence of proton ingestion. Lower panel: temporal 
evolution of the H-burning and He-burning contributions to luminosity. 



period. The final surface composition is therefore a 
combination of two different nucleosynthetic channels: 
the first being a consequence of the PIE and the second 
coming from the s-process in the radiative ^"^C pockets. 



3 Nucleosynthesis Results 

In Fig. [5] we report the evolution of selected heavy 
isotopes during the PIE: (representative of light. 
Is, s-process elements), ^^^La (representative of heavy, 
hs, s-process elements) and '^"^Pb (which is strongly 
produced by the s-process in metal-poor stars). Trac- 
ing s-process eleme nts have been selected basing on 
IBusso et al. I (j200 . During the first phases of proton 
ingestion (see panel a of Fig. O, the '^^C abundance 
(about 10~^ in mass fraction) does not allow an ap- 
preciable s-process production: heavy elements abun- 
dances are therefore the initial ones (in the previous 
phases of the stellar evolution no neutron sources have 
been in fact efficiently activated) . During the following 
phase, up to the splitting event, we obtain a non neg- 
ligible production of heavy elements, with particular 
efficiency in the production of hs elements (see panel 
b of Fig.lS}. Once the splitting has occurred, the nucle- 
osynthesis in the two shells follow completely different 
behaviors. In the lower shell, where ^"^C is destroyed 



at high temperatures {Tmax = 2.5 x 10* K), lead and 
hs elements are produced. In the upper shell, instead, 
'^^C burns at lower temperatures {Tmax < 10* K, see 
panel c and d in Fig[l} and it is continuously added 
by the active proton burning. The final effect is a net 
production of Is elements only (see panel d in Fig. [S}. 
When the envelope dredges up the isotopes synthe- 
sized in the upper shell, the surface distribution of the 
model is therefore enriched in Is elements, showing hs 
and lead abundances close to the initial ones (blue dot- 
ted curve in Fig. The mass dredged up during the 
d-TDU is AMtdu ~ 1.7 x 10"^ Mq; the correspond- 
ing efficiency (defined by means of the A factoiQ) is 
2.7. After the d-TDU episode, the model follows a 
standard evolution. The isotopes synthesized in the 
lower shell are diluted in the convective shell gener- 
ated by the following TP and their signature appears 
at the surface after the second TDU episode (red short- 
dashed curve in Fig. [S}. The surface s-process distri- 
bution after the 2""* TDU mimics the expect ed dis- 
tribution of models at low metallicity (Bistcrzo et al. I 
I2OO8), characterized by a growing sequence Is-hs-lead. 
Such a trend is even more strong as the number of 
TPs grows, displaying a final distribution character- 



is defined as the ratio between the mass of the H- 
depleted material that is dredged up and the mass of the 
material that has been burned by the H-shell during the 
previous interpulse period. 
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Figure 5: Temporal evolution of selected heavy isotopes during the PIE. 



ized by [ls/Fe]=2.02, [hs/Fe]=2.55 and [Pb/Fe]=3.50 
(dark solid curve in Fig. [B]). In the plotted final distri- 
bution we took into account the decays of radioactive 
isotopes (as, for example, the /3 decay of long-hved 
®^Zr to ^'^Nb) and, therefore, this distribution can be 
compared with observational data of extrinsic AGBs. 

3.1 Light elements 

As outlined in Section [21 a huge N enhancement re- 
sults as a consequence of proton ingestion from the 
top of th e convective shel l, a po ssibility early recog- 
nized bv iHoUowell et al. I (jl990h . After the d-TDU 
following the PIE, the surface is in fact enriched not 
only in ^'^C (as in a standard TDU), but also in ^''C 
and "N. In Fig. [7] we plot the [C/N] ratio as a func- 
tion of [C/Fe], while in Fig. [5] we report the loga- 
rithm of the ^^C/^'^C isotopic ratio as a function of 
[C/Fe]. In order to highlight the effects of the PIE, 
we report the same ratios for a mo del with initial 
mass M = 2 M0 and Z = 1 x 10"* ()Cristallo et al. I 
120091 ') ■ which follows a standard AGB evolution. More- 
over, we also display spectroscopic data of CEMP- 
s stars (both Giants and MS Stars) e xtracted from 
the SAGA database (|Suda et al. Il2008l 'l: our selecting 
criteria have been metallicity (-2.9<[Fe/H]<-2.2), a 
consistent s-elements enhancement ([Ba/Fe]>1.0) and 
data completeness (values of C, N and ^^C/^^C col- 
lected from single surveys). The final [C/N] ratio at- 
tained by our 1.5 M0 model better agrees with ob- 
servational data of CEMP-s stars, which show values 
ranging between and 1.3. For the sake of clarity, 
observational data show lower [C/Fe] ratios with re- 



spect to our model, but sinc e the majority of thes e 
stars belong to binary systems (|Lucatello et al. 1l2005l ) , 
a dilution factor has to possibly be accounted for (de- 
pending on the masses and the geometry of the system) 
(jAoki et al. ll2008l : lThompson et al. II2008D . Moreover, 
other physical mechanisms, such as gravitational set- 
tling, thermohaline mixing and the occurrence of the 
First Dredge Up, eventually led to a further dilution 
(Stancliffc & Glcbbcck (2008) and references therein) . 
Most of CEMP-s stars also show very low ^^C/^^C 
ratios (for some of them even close to the equilib- 
rium value), therefore proving that the material within 
their envelopes experienced some H-burning. Our final 
value (~ 80) is strongly reduced with respect to mod- 
els without the PIE (in which the ^^C/^^C ratio is 
larger than 10000) and, therefore, at low metallicities 
PIE help explaining the low observed carbon isotopic 
ratios. Observational data reported in Fig. [5] confirm 
this trend. Note that masses lower than 1.5 Mq expe- 
rience fewer TDU episodes and, therefore, lower final 
i2(-i^13q j^j.g expected with respect to the model pre- 
sented in this paper (Cristallo et al., in preparation). 
Moreover, we remember that large enough masses (the 
exact value depending on the meta Uicity) experienc e 
the Hot Bottom Burning (HBB) (jSueimoto I Il97ll ). 
showing therefore ^''C and ^'*N enriched surfaces. How- 
ever, the PIE can not occur at meta llicities larger than 
[Fe/H]> -2 (|Cristallo et al. II2007I ) and the mass ex- 
periencing HBB grows when increasing the metallicity. 
Additional extra-mixing mechanism, working during 
both the RGB and the AGB phase, are therefore re- 
quested to reproduce the low ^^C/^^ C ratios detected 
in stars of intermediate metallicity fsee lLebzelter et al. I 
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Figure 6: Surface abundances (in spectroscopic notation) of a model with initial mass M=1.5 M0 model 
and Z = 5 X 10~^ after selected TDU episodes. Gaps at Z = 43 and Z = 61 correspond to unstable 
elements, Tc and Pm, respectively. 



120081 1 as well in Galactic Open and Globular Clusters 
(see, e.g.. lGratton et al. l[2000l 'l. 

Finally, let us discuss a further interesting feature, 
already noted by Iwamoto et al. I (|2004). marking the 
PIE: a very large production of ^Li. We obtain a fi- 
nal log(£(Li)) = 3.74, about two order of magnitude 
highe r than the lithium pla teau observed at low metal- 
licity (|Spite &i Spite lll982t ). The ^Li production mech- 
anism is very simple: during proton ingestion, ^He nu- 
clei are ingested within the convective shell (note that 
the envelope presents a quite large '^He abundance, 
X('^He) ^ 3 X lO"*, due to the occurrence of the First 
Dredge Up). During the PIE, '^He is rapidly captured 
by the abundant *He, leading to the synthesis of ^Be. 
This isotope then captures an electron and synthesizes 
''Li. At high temperatures, lithium has a very large 
cross section against proton captures and, if created in 
an hot environment, is rapidly destroyed. However, 
given the rapidity of the stellar evolution after the 
splitting event, the freshly synthesized beryllium has 
not enough time to capture electrons in a hot environ- 
ment and, therefore, is rapidly diluted in the envelope 
by the d-TDU. 



ingestion has on the physical stellar structure and on 
the following nucleosynthesis. As a by-product of the 
PIE, we found a large production of ^'^C, ^"^N and ^Li. 
Moreover, we discussed the surface s-process distribu- 
tions after the deep TDU following the PIE and after 
the standard TDU episodes. We are now exploring 
the effects of the Proton Ingestion Episode on mod- 
els with different masses (0.85 < M/Mq < 2.0) and 
Z = 5 X 10~^ (Cristallo et al., in preparation). Pro- 
ton ingestion occurs in all these models. While in the 
1.5 and 2.0 Mq mass models the effects of proton in- 
gestion on the final surface compositions are integrated 
by the occurrence of standard TDU episodes (as shown 
in this contribution), lower masses (0.85 and 1.0 M0) 
show peculiar final surface compositions, determined 
by the PIE only. Note that PIE remarkably lowers the 
minimum mass for the occurrence of TDU. 
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[C/Fe] 



Figure 7: [C/N] ratio as a function of [C/Fe] for 
the M = 1.5 Mq model with Z = 5x IQ-^. We 
also report a model with M = 2.0 M0 and Z — 
1 X 10~^. As a comparison, we plot spectroscopic 
data of Giant and MS CEMP-s stars. See text for 
details. 




[C/Fe] 



Figure 8: ^^C/^^C isotopic ratio as a function of 
[C/Fe] for the M = 1.5 Mq model with Z = 5x 
10"^. We also report a model with M — 2.0 M© 
and Z — 1 X 10^^. As a comparison, we plot 
spectroscopic data of Giant and MS CEMP-s stars. 
See text for details. 
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